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Abstract 

In the context of the usual calibration model, we consider the case 
in which the independent variable is unobservable, but a pre-fixed 
value on its surrogate is available. Thus, considering controlled vari- 
ables and assuming that the measurement errors have equal variances 
we propose a new calibration model. Likelihood based methodology 
is used to estimate the model parameters and the Fisher information 
matrix is used to construct a confidence interval for the unknown value 
of the regressor variable. A simulation study is carried out to asses the 
effect of the measurement error on the estimation of the parameter of 
interest. This new approach is illustrated with an example. 
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1 Introduction 

In the first stage of a calibration problem, a pair of data sample (xj, Yi), i = 
1, 2, • • -n is observed. In the second stage, it is observed one or more values, 
which are the responses corresponding to a single unknown value of the re- 
gressor variable, X . The first and second stage equations of the usual linear 
calibration model are defined, respectively, as 

Yi = a + Pxi + ei, i = 1,2 ■••,n, (1.1) 
Y 0i = a + pX + e u i = n + l,n + 2,---,n + k. (1.2) 

It is considered the following assumptions: 

• xi, X2, ■ ■ ■ , x n take fixed values, which are considered as true values. 

• e 1; e 2 , • • • , e n+ k are independent and normally distributed with mean 
and variance of. 

The model parameters are a, (3, X and of and the main interest is to estimate 
the quantity X . 

The maximun likelihood estimators of the usual calibration model are 
given by 

& = Y-$x, $ = ^, X = ^^, (1.3) 

i n n+k 

^ = — r[E(*-«-fe) 2 + E(^-io) 2 ], (1.4) 

71 "+" K i=l i=n+l 

where 

in i n i n 

x = -5>i, F = -]Ty, s- xy = -2( a ; i -x)(y i -y), 
n f^i n f^i n f^i 

j n n+fc 

In [10] an approximate expression is derived for the variance of the esti- 
mator Xq, which is derived through the propagation error law. Another 
approximation for the variance of Xo is given by the Fisher information of 
6 = (a, /3,X ,<r 2 ) which, after some length algebraic manipulations, it can 
be shown to be given by 
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The maximum likelihood estimator of 6 = (a, (3, X Q , of ) has approxi- 
mately normal distribution with mean and covariance matrix I(8)~ l , when 
k = qn, q 6 Q + and n — > oo. Thus, the approximation of order n~ l for the 
variance of X is given by 



<7. 2 



1 , 1 + (X - X ) 



(1.6) 



On the other hand, in [4] the size k of the second stage is considered fixed, 
so that expanding X in Taylor series around the point (a,/3) and ignoring 
terms of order less than n~ 2 , we can find the following approximations for 
the bias and variance of X , respectively, 
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In order to construct a confidence interval for X , we consider that 

Xq — Xq d 



V(X ) 



iV(0,l), (1.9) 



where V(X ) is the estimated variance computed according to (1.6) or (1.8). 
Hence, the approximated confidence interval for X with a confidence level 
(1 — a), is given by 

(X - z°yJv{X ),X + z^V(X )^j , (1.10) 
where z« is the quantile of order (1 — f ) of the standard normal distribution. 

The usual calibration problem has been discussed in the literature for 
several decades (see [l]-[6]). An illustration of this model is presented for 
example in [7] . We can find a review of the literature on statistical calibration 
in [8], where some approaches to the solution of the calibration problem are 
summarized. 

This model encounters applications in different areas, but it is not well 
suited in some instances as, for example, in chemical analysis, where the 
preparation process of standard solutions are subject to measurement error 
([10]). 

There exists some situations, as mentioned above, where the independent 
variable, Xi, is measured with error. In this case, [11] defines two types of 
observations: controlled and uncontrolled. 

In the uncontrolled situation, the usual procedure to obtain the true value 
of the independent variable generates an error and the observed value is 

Xi=Xi + 5i, i = l,---,n. (1-11) 
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We have that CC^ IS clll unknown quantity, Si is a measurement error and 
Xi is a random variable. Assuming that X{ IS cL parameter the model defined 
by (1.1) and (1.11) is named as functional model ([12]). In this case there 
exists correlation between the model error and the variable Xi. Assuming 
that Xi IS db random variable the model (1.1) and (1.11) is called as structural 
model ([12]). On the other hand, the model defined by (1.1), (1-2) and (1.11) 
is called as the functional or structural calibration model if Xt is assumed as 
a parameter or a random variable, respectively ([13]). 

The controlled observation is defined by a pre-fixed value Xi according 
to the experimenter convenience and a procedure is established in order to 
attain the pre-fixed value. The experiment gives the unobserved Xi and it is 
such that 

Xi = Xi-5i, i = l,---,n. (1.12) 

In this case, the fixed quantity is Xi, the measurement error is 5i and xi 
is the random variable. The model (1.1) and (1.12) is known as Berkson 
regression model ([9]). Notice that the model error and the quantity Xi are 
independent. The model defined by (1.1), (1-2) and (1.12) has not been 
considered before in the measurement error literature and in this work it will 
be called as the controlled calibration model. 

In the calibration model defined by (1.1), (1.2) and (1.11), the values of 
the regressor, Xi, from the first stage are randomly generated, whereas in 
the controlled calibration model, (1.1), (1.2) and (1.12), they are assumed as 
pre-fixed by the experimenter. 

This work is organized as follows. In Section 2, we derive the maximum 
likelihood estimators of the homoscedastic controlled calibration model by 
considering both cases: <r| unknown and known. In Section 3, a simulation 
study is undertaken to investigate the sensitivity of parameter estimates of 
the proposed model. In Section 4, an example is presented to illustrate our 
new approach. In Section 5, the concluding remark is presented. 

2 Parameter estimation 

In this section we study the controlled calibration model. From the equations 
(1.1), (1.2) and (1.12) we can write 

Yi = a + pXi + iet- P6i), i = l,2---,n, (2.1) 
Y 0i = a + [3X + 6i, i = n + l,n + 2,---,n + k. (2.2) 

with the following assumptions for the random errors 

• €i are independent N(0, of) random variables. 
. E(^)=0, V(5i)= a\. 

• cov(^, Sj)=0 for any i ^ j. 
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• cov(ej, Sj)—0 for all 

Some comments are in order here. The variable in (2.1) is controlled 
and the error model (e, — (3 Si) is independent of X^. The error model in (2.2) 
is only in function of error measure e« related to Y 0i , this model assume that 
there is not error in the preparation sample related to parameter X . We 
define the homoscedastic controlled calibration model by considering that the 
errors Si are independent and normally distributed with mean and constant 
variance, a 2 . The study of this model is carried out following similar analysis 
to the usual calibration model as summarized above. 

The maximum likelihood estimator for the homoscedastic controlled cal- 
ibration model is derived in the following. The logarithm of the likelihood 
function is given by: 

l(a, (3, X , a 2 , aj) oc -^log(a 2 + L3 2 aj) - -log(a 2 ) - 



J2(Yi-a-pXi) 2 + 



i n+k 

- ( Y oi - a ~ P*o) 



e i=n+l 



(2.3) 



Solving dl/da = and dl/dX = we have the maximum likelihood esti- 
mator of a and X , which are given, respectively, by 



a = Y - PX and X 



Y -a 



(2.4) 



From (2.3) and (2.4), it follows that the likelihood for (f3, a 2 , a 2 ) can be 
written as 



1(0,^) « \log{ol + fo])-\lo 9 {o 2 )-\ 



i n+k 

P(x l -x)r+- e (y 0i 



(2.5) 



,2 ^ , v y o) 2 
'e i=n+l 

Next, we consider two cases for a 2 . Firstly, we obtain the maximum like- 
lihood estimator of /3, a 2 and a 2 from (2.5). In the second case we assume 

?! 

for (3 and a 2 . 



that the variance a 2 is known and obtain the maximum likelihood estimators 



Case 1: unknown variance erf 

Taking the partial derivative of (2.5) with respect to f3, of and a 2 and 
equating to zero we obtain, respectively, 



(3a 2 s (a 2 e + (3 2 &l - S YY + (3S XY ) = (S XY - (3S xx )a 2 e , 



(2.6) 



5 



a 2 + (3 2 a 2 = S YY -2[3S XY + (3 2 S XX , 



(2.7) 
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where S xx 
i 

n 



U=i(*i ~ X) 2 ,S XY = ±E? =1 (Xi - X)(Yi - Y),S YY 



\ E?=i(^ - Y? and S YoYo = \ E?=n+i(Yoi ~ Y ) 2 , and the relevant estimator 
notation has been introduced. From (2.6) and (2.7) we have the following 
equations: 



((3S XX - S XY )(S YY - 2(3S XY + f3 2 S xx ) = 0, 



hence 



(3S XX — S XY — 
S YY -2PS XY + P 2 S XX = 0. 



or 



(2.9) 
(2.10) 



Therefore, from (2.9), we have that (3 = S XY /S XX . But, according to the 
Cauchy-Schwarz inequality, S XX S YY > S XY , hence (2.10) has real roots if 
and only if Yi = cXi, where c is a constant. 

The estimator of a 2 can be obtained from the equation (2.7) 



YY 



2PS XY + (3 2 S XX ) 



(3 2 

Likewise, from equations (2.7) and (2.8) we obtain the estimator of the 
variance a 2 

a 2 = S YoYo . (2.11) 

In order to find the variance of Xo, we need to derive the Fisher infor- 
mation matrix of 9 = (a, (3, X , a 2 , of), which can be shown to be given 
by 
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where 



l = P 2 a 2 s +a 2 . (2.12) 

When k = qn, q G Q + and n — > oo, the estimator 9 is approximately 
normally distributed with mean 9 and variance I(9)~ l , thus we have that the 
approximate variance to order n" 1 for X is given by 



Vi(Xo) 



^1 



k 



+ 



rial 



+ 



7 (X - X 



n5 



(2.13) 
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Considering k fixed and expanding X in a Taylor series around (a,/3) 
and ignoring terms of order less than n~ 2 , it can be shown that the bias and 
variance of X (the proof is given in Appendix A), are given by 

7 (X - X ) 



Bias(X ) 



VAX{)) = f 2 
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(2.14) 
(2.15) 



nkf3 2 S xx _ 

We can observe that the estimator of X is biased, but it is asymptotically 
unbiased. 

With relation to the variance of the estimator X , let us notice that when 
k = qn, q G Q + , and ignoring the terms of order less than n^ 1 the variance 
in (2.15) coincide with the variance given in (2.13), which was found through 
the Fisher information. Equation (2.13) consider large sample sizes in the 
first and second stage (n and k), whereas (2.15) consider large sample sizes 
in the first stage and a fixed sample size in the second stage. 

Notice that when erf = 0, (2.13) and (2.15) coincide with (1.6) and (1.8) 
of the usual model, respectively. 

Caso 2: known variance ol 



Assuming now that erf is known and equating to zero the partial derivative 
of (2.5) with respect to the parameters (3 and erf, we have the following 
equations, respectively, 

Pa 2 (a 2 e + P 2 a 2 s - S YY + PS XY ) = (S XY - pS xx )a 2 and (2.16) 
kS YoYo k n S YY -2(JS XY + P 2 S XX 



a 2 + /3Vf 



(2.17) 



The estimates of (3 and err are obtained using some iterative method to solve 
(2.16) and (2.17). 

Similarly, as in Case 1, the Fisher information matrix of 9 = (a, (3, X , erf ) 
is given by 
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where 7 is defined in (2.12). 

The large sample variance of X follows by inverting the Fisher informa- 
tion matrix and is given by 

V(X ) 



P 2 



1 7 7 ^ 



(2.19) 
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where, 



nX 2 a\ + £X V - 2nX Xer £ 4 - 2kX X 1 2 + nX 2 er 4 + fcX 2 7 2 
(ner 4 + k'y 2 ) Er=i X 2 + 2nkp 2 ~iaj - n 2 X 2 af - nkX^ 2 ' 

Notice that if erf = 0, the expression (2.19) is reduced to (1.6). 

To construct a confidence interval for X , for both cases erf unknown 
and known, we consider the interval (1.10), where V(X c) is the estimated 
variance that follows from (2.13), (2.15) or (2.19). 

3 Simulation study 

In this section we present a simulation study for both cases of the ho- 
moscedastic controlled calibration model: erf known and unknown. The 
objetive of this section is to study the performance of the estimators of the 
proposed model (Proposed-M) and verify the impact by considering errati- 
cally the usual model (Usual-M). 

It was considered 5000 samples generated from the homoscedastic con- 
trolled calibration model. In all samples, the value of the parameters a and 
/3 were 0.1 and 2, respectively. The range of values for the controlled vari- 
able was [0,2]. The fixed values for the controlled variable were x\ = 0, Xi = 
Xi-i + ^-j-, i — 2, • • • , n, and the parameter values X were 0.01 (extreme in- 
ferior value), 0.8 (near to the central value) and 1.9 (extreme superior value). 
It was considered erf = 0.04 and the parameter values of erf were 0.01 and 
0.1, which are named, respectively, as small and large variances. For the 
first and second stages we consider the sample of sizes n — 5, 20, 100 and 
k = 2, 20, 100, respectively. 

The empirical mean bias is given by J2^=i(X — X )/5000 and the empir- 
ical mean squared error (MSE) is given by T,f=i(X -X ) 2 /5000. The mean 
estimated variance of X is given by Y,f=i ^(X )/5000, with V(X ) = V^Xq) 
or V^Xo), where Vi(X ) is the estimated variance of (1.6), (2.13) or (2.19) 
and V2(Xq) is the estimated variance of (2.15). The theoretical variances of 
Xo denoted as Vi(Xq) and V^Xo), are referred, respectively, to the expres- 
sions (1.6), (2.13) or (2.19) and (2.15) evaluated on the relevant parameter 
values. In Appendix B it is presented the simulation results. 

Tables Bl, B2, B5 and B6 present the empirical bias, the empirical mean 
squares error, the theoretical variance and the estimated variance of Xo. 
In these tables, it is considered only the variance (1.6) of the usual model, 
because based on a simulation study in [15] it was shown that the variances 
(1.6) and (1.8) give similar results. 

Tables B3, B4 and B7 present the covering percentages and the confidence 
interval amplitudes constructed with a 95% confidence level for the parameter 
X . In Table B3, the covering percentages %i and % 2 and amplitudes A\ and 
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A 2 are referred to the confidence intervals constructed using the equations 
(2.13) and (2.19). 

Tables B1-B4 consider the homoscedastic controlled calibration model 
assuming that erf is unknown. 

In Table Bl the empirical bias and MSE of Xo are little and an addition in 
the size of the variance erf, described in Table B2, causes an increasing in the 
bias and MSE. Moreover, we have that the bias and MSE of X are smaller 
when X is near to the center value of the variation interval of the variable 
X. These tables show that for all n, k and X , the theoretical variances 
obtained using the expressions (2.19) and (2.15) are equal. This fact occurs 
also for the mean estimated variances. We verify also that when n > 20 and 
k > 20 the theoretical variances and the mean estimated variances from the 
proposed model are approximately equal. Observing these tables, we can also 
notice that there exists differences between the mean estimated variances of 
the usual and proposed models. 

Analyzing Tables B3 and B4, we observe that for all n and X when it 
is adopted erratically the usual model, the amplitudes decrease very much 
as the size of k increases. This causes the covering percentage to decrease 
moving away from 95%. Whereas, adopting the proposed model it is observed 
that when k increases the confidence interval amplitude decreases, but the 
covering percentages increase approaching 95%. Notice that the covering 
percentage % 1 and % 2 and the amplitudes A 1 and A 2 are approximately 
equal, the amplitudes are very small for X = 0.8. In these tables, we observe 
that when k = 20 or 100 and when n increases the amplitudes of the intervals 
decrease and the covering percentages approaches 95%. In most cases, the 
covering percentage obtained through the proposed model are greater than 
that for the usual model results and are close to 95%. 

Tables B5 and B7 describe the results for the controlled homoscedastic 
calibration model with erf known. The iterative method Quasi-Newton [14] 
has been used. 

In Tables B5 and B6 we have that the empirical bias and SME decrease as 
the size of n or k increase and they are small when X Q is near to the central 
value of the variation interval, X = 0.8. When erf is small (Table B5), for all 
n and k, the empirical values of MSE from the usual and proposed model are 
close to the theoretical variance, but only the mean estimated variance from 
the proposed model is close to the theoretical variance. When erf is large 
(Table B6), in general, the empirical MSE and the mean estimated variance 
from the usual and proposed model are different, but the values supplied by 
the proposed model are very close to the theoretical variance. 

Analyzing Table B7, we can make similar comments to the ones we made 
about Tables B3 and B4. 
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4 Aplication 



In this section we test our model, considering both cases as known and 
unknown, using the data supplied by the chemical laboratory of the " Instituto 
de Pesquisas Tecnologicas (IPT)" - Brasil. We also consider the usual model 
in order to observe the performance of the proposed model. Our main interest 
is to estimate the unknown concentration value X of two samples A and B 
of the chemical elements cromo and cadmium. 

Tables 1 and 4 present the fixed values of concentration of the standard 
solutions and the corresponding intensities for the cromo and cadmium ele- 
ment, respectively, which are supplied by the plasma spectrometry method. 
This data is referred to as the first stage of the calibration model. 

Tables 2 and 5 present the intensities corresponding to 3 sample solutions 
from the sample A and B. This data is referred to as the second stage of the 
calibration model. 

Tables 3 and 6 describe the estimates of a, (3,X , V(X ), <rj and the 
confidence interval amplitude U(X ) from the homoscedastic controlled cal- 
ibration model of the samples A and B for the chemical elements cromo and 
cadmium. The values of the variance a\ considered as known are obtained 
from an external study carried out by the IPT, which are a\ = 2, 5865-E — 06 
for the cromo element and <rf = 0.0017E 1 + 02 for the cadmium element. As 
seen in Section 2, in order to obtain the estimates of the parameters (3 and 
<y\ of the proposed model when erf is known, iterative methods are required. 
In order to solve the system of equations (2.16) and (2.17) it was used the 
Quasi-Newton iterative method. It is also presented the estimates from the 
usual model. The estimates of the variance of X Q are computed using the 
relevant expressions (1.6), (2.13) or (2.19). The amplitude U(X ) is given by 
the product of the squared root of the estimated variance of X and 1.96. 

In Tables 3 and 6 we can observe that the estimates of a and (5 supplied 
by the usual model is equal to the proposed model when erf is unknown and 
they are equal for samples A and B, this occurs because the expression of the 
estimators a and j3 of both models are equal and they only depend on the first 
stage of the calibration model. These estimates are slightly different when 
compared with the estimates from the proposed model when cr| is known. 
With respect to the estimate of X , we observe that there is no difference of 
the estimates supplied by the usual and the proposed models of the cromo 
and cadmium element in both samples A and B, respectively. The estimates 
of the concentration of the sample A, of the elements cromo and cadmium, 
are outside of the variation range of the standard solution concentrations. We 
verify that, except to the sample B of the cadmium element, the estimates 
of the variance of X and the amplitude U(X ) from the usual model are 
greater than the estimates supplied by the both proposed models. 
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Table 1: Concentration (mg/g) and intensity of the standard solutions of 
cromo element. 



x l 


Intensity 


0,05 


6455,900 


0,11 


13042,933 


0,26 


32621,733 


0,79 


97364,500 


1,05 


129178,100 



Table 2: Intensity of the sample solutions A and B of cromo element. 



Intensity 


Sample A 


Sample B 


1465,0 
1351,0 
1495,6 


10173,6 
10516,9 
10352,2 



Table 3: Estimates of a, (3, X , V(X ) and the confidence interval amplitude 
U(X ) from the usual and proposed model for the samples A and B of cromo 



element. 



Parameters 


Sample A 


Sample B 


Usual-M 


Proposed-M 


Usual-M 


Proposed-M 


unknown of 


known of 


unknown of 


known of 


a 


123,574 


123,574 


123,889 


123,574 


123,574 


124,021 




l,23E+05 


l,23E+05 


l,23E+05 


l,23E+05 


1.23E+05 


l,23E+05 


X 


0,011 


0,011 


0,011 


0,083 


0,083 


0,083 


V(X ) 


9.80E-07 


9.15E-07 


1.35E-06 


l,16E-06 


1.13E-06 


1.71E-06 


-I 




1.60E-06 






5.48E-07 




U(X ) 


2.55E-03 


2.46E-03 


2.99E-03 


2J7E-03 


2.73E-03 


3.36E-03 



Table 4: Concentration (mg/g) and intensity of the standard solutions of 
cadmium element. 



x l 


Intensity 


0,05 


4,89733 


0,10 


9,706 


0,25 


23,41333 


0,73 


69,73 


1,01 


96,85667 
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Table 5: Intensity of the sample solutions A and B of cadmium element. 



Intensity 


Sample A 


Sample B 


0,679 
0,6837 
0,6846 


5,066 
5,027 
5,085 



Table 6: Estimates of a, f3, X , V(X ) and the confidence interval amplitude 
U(X ) from the usual and homoscedastic models for the samples A and B of 
cadmium element. 



Parameters 


Sample A 


Sample B 


Usual-M 


Proposed-M 


Usual-M 


Proposcd-M 


unknown <r^ 


known crjj 


unknown 


known a£ 


a 


-0,156 


-0,156 


-0,158 


-0,156 


-0,156 


-0,158 


(1 


95,828 


95,828 


95,831 


95,828 


95,828 


95,831 


x 


8J5E-03 


8.75E-03 


8.77E-03 


0,054 


0,054 


0,054 


V(X ) 


4,06E-06 


3.72E-06 


1.26E-06 


3.81E-06 


3,32E-06 


1.17E-06 






8.31E-06 






8,24E-06 




U(X ) 


5,18E-03 


4.96E-03 


2.89E-03 


5.02E-03 


4,68E-03 


2.78E-03 



5 Concluding remarks 

In general, the simulation study reveals that the proposed model is sensible 
to the presence of error related to the independent variable and gives better 
results in contrast to the usual model results. It was noticed that when the 
error variance erf increases, the mean estimated variance of X obtained us- 
ing the usual model moves away from the theoretical value. In the example 
above, the confidence interval amplitude from the proposed models are sup- 
plied by the incorporation of error due to the lecture of equipment and the 
preparation of the standard solutions. It is observed that despite the classical 
model only considers the error originated from the lecture of the equipment, 
the amplitude is greater than the obtained by the new approach. 
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Appendix 



A Bias and variance for the maximum like- 
lihood estimator 

In the following we derive the bias (2.14) and the variance (2.15) of the 
estimator Xo from the homoscedastic controlled calibration model when erf 
is known. 

Considering the model (2.1) and (2.2), the estimator X = (Y — a)//3 
can be expressed as 

4=J+ w I W i (A1) 

ft 

where e = E?= + „+i and = E? = i(q - pSJ/n. 

Considering k fixed, expanding 1/(3 in a Taylor series around (3 and ig- 
noring terms of order less than n~ 2 , we obtain the expected value of (A.l), 
given by 

£ ™ = A '° + 1 §^ 21 - < A - 2 > 

From this last equation we get the bias (2.14). 

To derive the variance (2.15) we take the variance of (A.l), which is given 

by 

V(X ) = (3 2 (X - XfV{l) + V{%) + V(t). (A.3) 

We call attention to the fact that (A.3) is only expressed as a function of the 
related variances because the corresponding covariances are zero. The vari- 
ances V(l/(3),V(e /(3) and V{<f>/(3) can be obtained by expanding 1//3, e //3 
and <f)//3 in a Taylor series around (3 and ignoring terms of order less than 
n~ 2 . They are given by 

V(l/P) = (A.4) 

V(e /P) = ^ + 3^V0), (A.5) 
Vm = ^- (A.6) 

Substituing (A.4), (A.5) and (A.6) in (A.3), then, the variance (2.15) is 
obtained. 
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B Tables 

Table Bl. Empirical bias and mean squared error, theoretical variance and 
the mean estimated variance of X , for erf = 0, 01 and unknown. 









Empirical 


Theoretical 


Mean of V(X ) 


Xq 


n 


k 






Proposcd-M 


Usual-M 


Proposed-M 








Bias 


MSE 


Vi(X ) 


V 2 (X ) 


Vi(Xo) 


Vi(Xo) 


V 2 {X ) 


0,01 


5 


2 


-0,0060 


0,0180 


0,0170 


0,0170 


0,0120 


0,0100 


0,0100 






20 


-0,0087 


0,0130 


0,0120 


0,0120 


0,0072 


0,0120 


0,0120 






100 


-0,0060 


0,0130 


0,0120 


0,0120 


0,0065 


0,0120 


0,0120 




20 


2 


-0,0038 


0,0086 


0,0087 


0,0087 


0,0120 


0,0052 


0,0053 






20 


-0,0028 


0,0043 


0,0042 


0,0042 


0,0033 


0,0040 


0,0040 






100 


-0,0032 


0,0038 


0,0038 


0,0038 


0,0022 


0,0036 


0,0036 




100 


2 


-0,0023 


0,0058 


0,0058 


0,0058 


0,0100 


0,0027 


0,0027 






20 


-0,0002 


0,0013 


0,0013 


0,0013 


0,0016 


0,0012 


0,0012 






100 


-0,0007 


0,0008 


0,0009 


0,0009 


0,0007 


0,0009 


0,0009 


0,8 


5 


2 


-0,0011 


0,0094 


0,0093 


0,0094 


0,0079 


0,0045 


0,0046 






20 


-0,0034 


0,0050 


0,0048 


0,0048 


0,0029 


0,0045 


0,0046 






100 


-0,0007 


0,0047 


0,0044 


0,0044 


0,0024 


0,0045 


0,0045 




20 


2 


0,0005 


0,0063 


0,0061 


0,0061 


0,0095 


0,0028 


0,0029 






20 


0,0007 


0,0016 


0,0016 


0,0016 


0,0015 


0,0015 


0,0015 






100 


-0,0001 


0,0012 


0,0012 


0,0012 


0,0008 


0,0012 


0,0012 




100 


2 


0,0005 


0,0050 


0,0052 


0,0052 


0,0099 


0,0021 


0,0021 






20 


-0,0001 


0,0007 


0,0007 


0,0007 


0,0011 


0,0007 


0,0007 






100 


-0,0003 


0,0003 


0,0003 


0,0003 


0,0003 


0,0003 


0,0003 


1,9 


5 


2 


0,0041 


0,0160 


0,0150 


0,0160 


0,0120 


0,0093 


0,0094 






20 


0,0026 


0,0110 


0,0110 


0,0110 


0,0065 


0,0110 


0,0110 






100 


0,0076 


0,0110 


0,0110 


0,0110 


0,0058 


0,0110 


0,0110 




20 


2 


0,0006 


0,0079 


0,0082 


0,0082 


0,0110 


0,0049 


0,0049 






20 


0,0040 


0,0039 


0,0037 


0,0037 


0,0030 


0,0035 


0,0035 






100 


0,0008 


0,0033 


0,0033 


0,0033 


0,0019 


0,0031 


0,0031 




100 


2 


0,0020 


0,0057 


0,0057 


0,0057 


0,0100 


0,0025 


0,0025 






20 


0,0003 


0,0012 


0,0012 


0,0012 


0,0015 


0,0011 


0,0011 






100 


0,0003 


0,0008 


0,0008 


0,0008 


0,0006 


0,0008 


0,0008 
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Table B2. Empirical bias and mean squared error, theoretical variance 
and the mean estimated variance of X , for a\ — 0, 1 and unknown. 









Empirical 


Theoretical 


Mean of V(X ) 




n 


k 






Proposcd-M 


Usual-M 


Proposed-M 








Bias 


MSE 


Vi(X ) 


V 2 (X Q ) 


Vi(X ) 


Vi(X ) 


V 2 {X ) 


0,01 


5 


2 


-0,0510 


0,1000 


0,0700 


0,0710 


0,0770 


0,0680 


0,0690 






20 


-0,0500 


0,0950 


0,0660 


0,0660 


0,0220 


0,0660 


0,0660 






100 


-0,0510 


0,0950 


0,0650 


0,0650 


0,0130 


0,0730 


0,0730 




20 


2 


-0,0180 


0,0280 


0,0250 


0,0250 


0,0670 


0,0240 


0,0240 






20 


-0,0160 


0,0230 


0,0210 


0,0210 


0,0140 


0,0210 


0,0210 






100 


-0,0170 


0,0230 


0,0200 


0,0200 


0,0055 


0,0210 


0,0210 




100 


2 


-0,0046 


0,0094 


0,0093 


0,0093 


0,0580 


0,0069 


0,0069 






20 


-0,0026 


0,0048 


0,0048 


0,0048 


0,0082 


0,0048 


0,0048 






100 


-0,0033 


0,0043 


0,0044 


0,0044 


0,0029 


0,0044 


0,0044 


0,8 


5 


2 


-0,0084 


0,0370 


0,0290 


0,0290 


0,0460 


0,0250 


0,0250 






20 


-0,0095 


0,0270 


0,0240 


0,0240 


0,0071 


0,0200 


0,0200 






100 


-0,0072 


0,0290 


0,0240 


0,0240 


0,0037 


0,0200 


0,0200 




20 


2 


-0,0030 


0,0120 


0,0110 


0,0110 


0,0530 


0,0085 


0,0086 






20 


-0,0040 


0,0068 


0,0066 


0,0066 


0,0061 


0,0064 


0,0064 






100 


-0,0031 


0,0063 


0,0062 


0,0062 


0,0017 


0,0060 


0,0060 




100 


2 


-0,0011 


0,0063 


0,0062 


0,0063 


0,0550 


0,0038 


0,0038 






20 


-0,0015 


0,0017 


0,0017 


0,0017 


0,0057 


0,0017 


0,0017 






100 


-0,0011 


0,0014 


0,0013 


0,0013 


0,0013 


0,0013 


0,0013 


1,9 


5 


2 


0,0430 


0,1090 


0,0630 


0,0630 


0,0830 


0,0750 


0,0750 






20 


0,0450 


0,0860 


0,0580 


0,0580 


0,0210 


0,0650 


0,0650 






100 


0,0410 


0,0860 


0,0580 


0,0580 


0,0110 


0,0600 


0,0600 




20 


2 


0,0160 


0,0260 


0,0230 


0,0230 


0,0650 


0,0210 


0,0210 






20 


0,0140 


0,0190 


0,0180 


0,0180 


0,0130 


0,0180 


0,0180 






100 


0,0170 


0,0200 


0,0180 


0,0180 


0,0048 


0,0180 


0,0180 




100 


2 


0,0050 


0,0088 


0,0087 


0,0088 


0,0570 


0,0063 


0,0063 






20 


0,0030 


0,0043 


0,0042 


0,0042 


0,0078 


0,0042 


0,0042 






100 


0,0020 


0,0039 


0,0038 


0,0038 


0,0026 


0,0038 


0,0038 
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Table B3. Covering percentage (%) and amplitude (A) of the intervals 
with a 95% confidence level for the parameter X , when a\ = 0, 01 and 
unknown. 



x 


n 


k 


Usual-M 


Proposcd-M 








% 


A 


%l 


Ai 


%2 


A 2 


0,01 


5 


2 


83,04 


0,40 


79,34 


0,36 


79,37 


0,36 






20 


84,95 


0,32 


91,15 


0,41 


91,15 


0,41 






100 


83,72 


0,31 


92,24 


0,42 


92,24 


0,42 




20 


2 


96,46 


0,42 


83,48 


0,28 


83,52 


0,28 






20 


90,19 


0,22 


92,47 


0,24 


92,47 


0,24 






100 


86,16 


0,18 


92,78 


0,23 


92,78 


0,23 




100 


2 


98,84 


0,40 


73,16 


0,19 


73,16 


0,19 






20 


96,71 


0,16 


94,14 


0,14 


94,14 


0,14 






100 


91,90 


0,11 


94,68 


0,12 


94,68 


0,12 


0,8 


5 


2 


85,43 


0,32 


74,33 


0,24 


74,39 


0,24 






20 


85,16 


0,21 


91,34 


0,26 


91,34 


0,26 






100 


85,04 


0,19 


92,50 


0,25 


92,50 


0,25 




20 


2 


97,90 


0,38 


73,55 


0,20 


73,55 


0,20 






20 


93,55 


0,15 


93,89 


0,15 


93,89 


0,15 






100 


86,53 


0,11 


93,12 


0,13 


93,12 


0,13 




100 


2 


99,41 


0,39 


65,05 


0,16 


65,05 


0,16 






20 


98,54 


0,13 


94,05 


0,10 


94,05 


0,10 






100 


94,56 


0,07 


94,86 


0,07 


94,86 


0,07 


1,9 


5 


2 


82,47 


0,39 


78,05 


0,35 


78,13 


0,35 






20 


84,50 


0,31 


90,92 


0,39 


90,92 


0,39 






100 


84,75 


0,29 


92,83 


0,39 


92,83 


0,39 




20 


2 


96,79 


0,41 


83,09 


0,26 


83,11 


0,26 






20 


91,32 


0,21 


93,29 


0,23 


93,31 


0,23 






100 


86,43 


0,17 


93,32 


0,22 


93,32 


0,22 




100 


2 


98,88 


0,40 


73,06 


0,19 


73,06 


0,19 






20 


97,12 


0,15 


94,31 


0,13 


94,31 


0,13 






100 


92,56 


0,10 


94,74 


0,11 


94,74 


0,11 
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Table B4. Covering percentage (%) and amplitude (A) of the intervals 
with a 95% confidence level for the parameter X , when erf = 0, 1 and un- 
known. 



X 


n 


k 


Usual-M 


Proposed-M 








% 


A 


%i 


Ai 


%2 


A 2 


0,01 


5 


2 


84,89 


0,94 


80,73 


0,85 


80,79 


0,86 






20 


64,10 


0,51 


82,75 


0,86 


82,75 


0,86 






100 


52,09 


0,39 


82,42 


0,86 


82,42 


0,86 




20 


2 


99,40 


0,99 


91,04 


0,58 


91,10 


0,58 






20 


87,10 


0,45 


92,60 


0,55 


92,62 


0,55 






100 


65,70 


0,28 


92,16 


0,55 


92,18 


0,55 




100 


2 


100,00 


0,94 


87,60 


0,32 


87,66 


0,32 






20 


98,82 


0,36 


94,88 


0,27 


94,88 


0,27 






100 


89,12 


0,21 


94,84 


0,26 


94,84 


0,26 


0,8 


5 


2 


90,27 


0,73 


80,63 


0,52 


80,75 


0,52 






20 


64,39 


0,31 


82,68 


0,51 


82,74 


0,51 






100 


50,57 


0,23 


84,08 


0,50 


84,08 


0,50 




20 


2 


99,92 


0,88 


87,30 


0,35 


87,50 


0,35 






20 


91,80 


0,30 


92,78 


0,31 


92,82 


0,31 






100 


67,46 


0,16 


92,38 


0,30 


92,38 


0,30 




100 


2 


100,00 


0,91 


77,88 


0,22 


77,98 


0,22 






20 


99,94 


0,29 


94,98 


0,16 


95,04 


0,16 






100 


94,60 


0,14 


94,92 


0,14 


94,92 


0,14 


1,9 


5 


2 


85,63 


0,91 


81,42 


0,81 


81,42 


0,81 






20 


61,71 


0,49 


81,89 


0,82 


81,89 


0,82 






100 


50,71 


0,37 


81,24 


0,81 


81,24 


0,81 




20 


2 


99,54 


0,97 


91,12 


0,54 


91,18 


0,55 






20 


88,10 


0,43 


92,74 


0,52 


92,76 


0,52 






100 


66,44 


0,26 


92,86 


0,51 


92,86 


0,51 




100 


2 


100,00 


0,93 


86,38 


0,30 


86,38 


0,30 






20 


99,12 


0,35 


94,76 


0,25 


94,76 


0,25 






100 


89,40 


0,20 


95,00 


0,24 


95,00 


0,24 
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Table B5. Empirical bias and mean squared error, theoretical variance 
and the mean estimated variance of X , for <rf = 0, 01 and known. 









Empirical 


Theoretical 


Mean of V(X ) 




n 


k 


Usual-M 


Proposcd-M 


Proposcd-M 


Usual-M 


Proposcd-M 








Bias 


MSE 


Bias 


MSE 


V(X ) 






0,01 


5 


2 


-0,0290 


0,0210 


-0,0280 


0,0210 


0,0170 


0,0180 


0,0160 






20 


-0,0290 


0,0140 


-0,0320 


0,0140 


0,0120 


0,0081 


0,0130 






100 


-0,0240 


0,0140 


-0,0270 


0,0140 


0,0120 


0,0070 


0,0130 




20 


2 


-0,0081 


0,0091 


-0,0064 


0,0090 


0,0086 


0,0130 


0,0076 






20 


-0,0060 


0,0043 


-0,0072 


0,0043 


0,0041 


0,0034 


0,0041 






100 


-0,0038 


0,0038 


-0,0060 


0,0038 


0,0037 


0,0022 


0,0038 




100 


2 


-0,0011 


0,0056 


-0,0005 


0,0056 


0,0058 


0,0100 


0,0053 






20 


-0,0002 


0,0013 


-0,0001 


0,0013 


0,0013 


0,0016 


0,0012 






100 


-0,0009 


0,0009 


-0,0012 


0,0009 


0,0009 


0,0007 


0,0009 


0,8 


5 


2 


-0,0074 


0,0110 


-0,0072 


0,0100 


0,0093 


0,0120 


0,0085 






20 


-0,0046 


0,0051 


-0,0051 


0,0052 


0,0048 


0,0032 


0,0049 






100 


-0,0076 


0,0048 


-0,0082 


0,0048 


0,0044 


0,0025 


0,0046 




20 


2 


-0,0034 


0,0063 


-0,0031 


0,0063 


0,0061 


0,0100 


0,0052 






20 


0,0001 


0,0016 


0,0000 


0,0016 


0,0016 


0,0015 


0,0016 






100 


-0,0009 


0,0012 


-0,0013 


0,0012 


0,0012 


0,0008 


0,0012 




100 


2 


0,0000 


0,0053 


0,0002 


0,0053 


0,0052 


0,0099 


0,0048 






20 


0,0000 


0,0007 


0,0000 


0,0007 


0,0007 


0,0011 


0,0007 






100 


-0,0001 


0,0003 


-0,0002 


0,0003 


0,0003 


0,0003 


0,0003 


1,9 


5 


2 


0,0200 


0,0180 


0,0200 


0,0180 


0,0150 


0,0170 


0,0140 






20 


0,0240 


0,0120 


0,0260 


0,0130 


0,0110 


0,0071 


0,0120 






100 


0,0200 


0,0130 


0,0230 


0,0130 


0,0100 


0,0062 


0,0110 




20 


2 


0,0037 


0,0082 


0,0021 


0,0081 


0,0082 


0,0120 


0,0071 






20 


0,0059 


0,0037 


0,0066 


0,0037 


0,0037 


0,0030 


0,0036 






100 


0,0033 


0,0032 


0,0051 


0,0032 


0,0032 


0,0020 


0,0033 




100 


2 


0,0020 


0,0058 


0,0015 


0,0057 


0,0057 


0,0100 


0,0053 






20 


0,0003 


0,0012 


0,0002 


0,0012 


0,0012 


0,0015 


0,0011 






100 


0,0003 


0,0008 


0,0006 


0,0008 


0,0008 


0,0006 


0,0008 
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Table B6. Empirical bias and mean squared error, theoretical variance 
and the mean estimated variance of X , for a\ = 0, 1 and known. 









Empirical 


Theoretical 


Mean of V(X ) 




n 


k 


Usual-M 


Proposcd-M 


Proposcd-M 


Usual-M 


Proposcd-M 








Bias 


MSE 


Bias 


MSE 


V(X ) 






0,01 


5 


2 


-0,4330 


0,5590 


-0,3830 


0,4580 


0,0590 


0,3650 


0,2310 






20 


-0,1140 


0,1310 


0,4140 


1,0880 


0,0540 


0,0250 


0,0760 






100 


-0,1890 


0,1540 


-0,0290 


0,7730 


0,0540 


0,0200 


0,1250 




20 


2 


-0,0930 


0,0430 


-0,0770 


0,0360 


0,0210 


0,0950 


0,0380 






20 


-0,0490 


0,0210 


-0,0510 


0,0210 


0,0160 


0,0150 


0,0180 






100 


-0,0430 


0,0190 


-0,0250 


0,0640 


0,0150 


0,0058 


0,0170 




100 


2 


-0,0200 


0,0110 


-0,0160 


0,0097 


0,0084 


0,0630 


0,0110 






20 


-0,0077 


0,0041 


-0,0085 


0,0039 


0,0037 


0,0084 


0,0038 






100 


-0,0066 


0,0037 


-0,0087 


0,0037 


0,0033 


0,0029 


0,0033 


0,8 


5 


2 


-0,0500 


0,0620 


-0,0430 


0,0550 


0,0280 


0,1150 


0,0620 






20 


-0,0450 


0,0380 


0,0820 


0,0810 


0,0240 


0,0086 


0,0320 






100 


-0,0530 


0,0510 


-0,0069 


0,0980 


0,0230 


0,0055 


0,0360 




20 


2 


-0,0280 


0,0150 


-0,0240 


0,0140 


0,0110 


0,0740 


0,0210 






20 


-0,0079 


0,0068 


-0,0085 


0,0069 


0,0064 


0,0065 


0,0067 






100 


-0,0055 


0,0065 


-0,0030 


0,0079 


0,0060 


0,0018 


0,0062 




100 


2 


-0,0030 


0,0067 


-0,0021 


0,0065 


0,0062 


0,0600 


0,0091 






20 


-0,0018 


0,0018 


-0,0018 


0,0017 


0,0017 


0,0058 


0,0017 






100 


-0,0021 


0,0013 


-0,0025 


0,0013 


0,0013 


0,0013 


0,0013 


1,9 


5 


2 


0,3480 


0,3770 


0,3070 


0,3000 


0,0530 


0,2850 


0,1760 






20 


0,1400 


0,1630 


-0,3310 


0,9370 


0,0490 


0,0300 


0,0780 






100 


0,0410 


0,0180 


0,0270 


0,0510 


0,0140 


0,0051 


0,0150 




20 


2 


0,0970 


0,0430 


0,0790 


0,0350 


0,0190 


0,0930 


0,0350 






20 


0,1400 


0,1630 


-0,3310 


0,9370 


0,0490 


0,0300 


0,0780 






100 


0,1670 


0,1300 


-0,0005 


0,7090 


0,0480 


0,0160 


0,1120 




100 


2 


0,0200 


0,0093 


0,0160 


0,0087 


0,0080 


0,0620 


0,0110 






20 


0,0120 


0,0039 


0,0130 


0,0037 


0,0033 


0,0080 


0,0035 






100 


0,0072 


0,0031 


0,0089 


0,0031 


0,0029 


0,0027 


0,0030 
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Table B7. Covering percentage (%) and amplitude (A) of the inter 
with a 95% confidence level for the parameter X , when <r| = 0, 01. and 
and known. 













0,01 








0, 1 






n 


k 


Usual-M 


Proposcd-M 


Usual-M 


Proposed-M 








% 


A 


% 


A 


% 


A 


% 


A 


0,01 


5 


2 


92,10 


0,51 


91,20 


0,48 


95,06 


1,89 


92,40 


1,49 






20 


87,32 


0,34 


95,18 


0,44 


63,38 


0,59 


64,47 


1,09 






100 


84,50 


0,32 


95,19 


0,44 


52,69 


0,48 


89,24 


1,18 




20 


2 


97,46 


0,43 


90,12 


0,33 


99,87 


1,17 


96,27 


0,71 






20 


91,00 


0,23 


94,03 


0,25 


92,76 


0,48 


94,48 


0,52 






100 


86,73 


0,19 


95,21 


0,24 


70,38 


0,30 


90,35 


0,51 




100 


2 


97,01 


0,43 


90,00 


0,33 


100,00 


0,98 


93,41 


0,40 






20 


92,34 


0,23 


95,26 


0,25 


99,77 


0,36 


95,09 


0,24 






100 


85,93 


0,19 


94,55 


0,24 


92,70 


0,21 


94,42 


0,23 


0,8 


5 


2 


94,33 


0,41 


89,28 


0,35 


98,77 


1,31 


97,54 


0,98 






20 


86,84 


0,22 


95,03 


0,27 


61,85 


0,34 


81,25 


0,69 






100 


85,65 


0,19 


95,56 


0,27 


49,29 


0,27 


87,82 


0,69 




20 


2 


98,35 


0,39 


88,44 


0,27 


100,00 


1,03 


94,09 


0,53 






20 


93,84 


0,15 


94,41 


0,15 


94,44 


0,32 


95,11 


0,32 






100 


85,98 


0,11 


94,32 


0,14 


70,38 


0,17 


92,41 


0,31 




100 


2 


98,08 


0,39 


88,22 


0,27 


100,00 


0,95 


92,25 


0,34 






20 


92,75 


0,15 


93,96 


0,15 


99,89 


0,30 


94,89 


0,16 






100 


87,17 


0,11 


94,60 


0,14 


95,00 


0,14 


94,88 


0,14 


1,9 


5 


2 


92,30 


0,50 


90,61 


0,46 


96,26 


1,78 


92,14 


1,45 






20 


86,96 


0,33 


95,13 


0,42 


63,58 


0,56 


67,24 


1,03 






100 


85,84 


0,30 


94,73 


0,42 


49,30 


0,46 


83,75 


1,21 




20 


2 


97,61 


0,43 


89,92 


0,32 


100,00 


1,17 


96,27 


0,70 






20 


91,60 


0,21 


94,27 


0,23 


94,26 


0,46 


94,95 


0,49 






100 


86,16 


0,17 


94,99 


0,23 


72,44 


0,28 


94,18 


0,48 




100 


2 


97,04 


0,43 


89,70 


0,32 


100,00 


0,97 


94,27 


0,39 






20 


91,21 


0,21 


93,85 


0,23 


99,59 


0,35 


95,18 


0,23 






100 


87,04 


0,17 


95,11 


0,23 


93,38 


0,20 


94,98 


0,21 
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